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To  go  with  the  HYCOM-based  global  ocean  prediction 
system  scheduled  for  transition  to  NAVOCEANO  in  FY06 
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Nesting  to  higher  resolution  coastal  domains  using 
HYCOM  and/or  NCOM 

•  I  mplies  a  vertical  remapping  from  global  HYCOM  to  the 
target  vertical  structure  of  the  coastal  domain 

•  The  nesting  scheme  must  accurately  represent  flow  regimes 
with  widely  different  dynamics  and  time  scales  (e.g.  in  shallow 
water,  over  the  continental  slope  and  in  deep  water,  with  all 
three  present  in  many  cases). 


Mississippi  Bight 


Current  Status  of  Nesting 
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Currently  off-line 

Boundary  info  comes  from  archive  files 
Exact  boundary  condition  for  depth 
averaged  component 
Relaxation  in  buffer  zone  forT,S,P,u,v 
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•  Updating  frequency  limited  by  archive  file  frequency 

•  Don't  need  to  know  nested  area  in  advance 


Off-line: 

•  Boundary  information  comes  from 
archive  files 

•  Updating  frequency  limited  by 
archive  file  frequency 

•  Don't  need  to  know  nest  area  in  advance 


On-line  (not  yet  implemented): 

•  Local  model  runs  simultaneously 
with  regional/ basin-scale  model 

•  Coupled  via  a  vertical  remapper 

•  Need  to  know  nested  area  in 
advance 
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Nested  SSH,  u,v  after  180  days  (0.1-9  day  e-folding;  T-S-p,u,v  relax) 
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sea  surf,  height  year  6.SQ  (jul  16)  [02.4H] 


fa  ye  r =01  u-veloc.  year  6.50  (jul  16)  [02.4H] 


layer=Q1  v-veloc.  year  6.50  (jul  16)  [02.4 H] 


layer=01  u-veloc.  year  6.50  (jul  16)  [01. 2H] 
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0.08°  Florida  Current  Nested  Region 


Original  SSH,  u,v  after  180  days 


0.08°  Florida  Current  Nested  Region 


0.08°  Florida  Current  Nested  Region 


6.02  (jan  22)  [01. 3H] 


Nested  section  after  6  days:  Nested  section  after  6  days: 

0.1-day  e-folding  T-S-p  relax  0.1-day  e-folding  T-S-p-u-v  relax 


u-velocity  merid.sec.  83.04w  year  6.02  (jan  22)  [01. 2H] 


temperature  merid.sec.  83.04w  year  6.02  (jan  22)  [01 .2H] 


i-velocity  merid.sec.  S3.Q4w  year  6.02  (jan  22)  [01. 3H] 
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0.08°  Gulf  of  Mexico  nested 
inside  of  0.08°  North  Atlantic 


0,08°  North  AcJa ndc 


3  months  after  restart 
Relaxation  to  T,S,p 
1-10  day  e-folding  time 
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Plan  to  due  3x  nest  or  —2.7  km 


Relaxation  buffer  zones 


Gulf  of  California  Nesting 

0.08°  Pacific  Ocean  Model  0.08°  Nested  Gulf  of  California 


1-10  day  e-folding  time,  relaxation  to  T,S,p,u,v 
21  days  after  initialization  Courtesy:  J  oe  Metzger 
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1/25°  (3.5  km)  East  Asian  Seas  HYCOM 
Nested  inside  1/6  °  North  Pacific  HYCOM 


open  boundary 
conditions  from  1/6° 
North  Pacific  HYCOM 


Nested  model  has  same  vertical 
structure  as  Pacific  Ocean  model 
(20  layers) 


•  Currently  off-line 

•  Boundary  info  comes  from  archive  files 

•  Exact  boundary  condition  for  depth 
averaged  component 

•  Relaxation  in  buffer  zone  for  T,S,P,u,v 


velocity  merid.sec.1 24. 52e  year  7.79  (o cl  14)  [01 .1 H] 


temperature  merid.sec.1  24.52e  year  7.79  (oct  1 4)  [01 .1  H] 


Yellow  Sea  flow 
reversal  with  depth 


Isopycnals  over 
shelf  region 


1/25°  HYCOM  East  Asian  Seas  Model 
(nested  inside  1/6°  North  Pacific  Model 

North-south  cross-section  along  124.5°E 


East  China  Sea 


blue=westward  flow 
red=eastward  flow 


density  front 
associated  with 
sharp  topo  feature 
(can’t  resolve  with 
sigma  coordinates) 


Snapshot  on  Oct.  14 


1/25°  HYCOM  East  Asian  Seas  Model 
(nested  inside  1/6°  North  Pacific  Model 


red=eastward  flow 
blue=westward  flow 


North-south  cross-section  along  124.5°E 


velocity  merid.sec.124.52e  year  8.28  (apr  12)  [01 .1 H] 


Vel  ocity 


Temperature 


z-levels  and  sigmas 
over  shelf  and  in 
mixed  layer 


Snapshot  on  Apr.  12 


Future  Plans 


Coastal  Ocean 


Slope  to  Shelf  Energetics 
And  Exchange  Dynamics 

(Jacobs,  Teague,  Hogan,  Arnone) 


Nesting  Studies 

(Hogan,  Kindle,  Wallcraft) 

• Develop  HYCOM  coastal  capabilities 


IAS  d  0.08 


Evaluate  coastal  HYCOM  and  NCOM 


Measuring  and  modeling  of 
processes  that  affect  cross¬ 
slope  exchanges  in  the  GoM 


Evaluate  coupling  and  boundary 
conditions 


Mississippi  Bight  Domain 


•  Designed  for  coastal  ocean  Modeling 

•  Hybrid  sigma-Z  vertical  coordinate 
(sigma  is  terrain- following) 

•  Transition  between  coordinates 
fixed  (ideally  at  shelf  break) 

•  Limited  representation  of  topog. 

(full  cells  in  z-level  mode) 

•  ~  3  times  faster  per  layer  than 
HYCOM 

•  Developed  from  POM  by  Paul  Martin 
(NRL) 


•  Designed  for  deep  water  domains 
and  accurate  deep  to  shallow  water 
transition 

•  Hybrid  isopycnal-sigma-Z 

•  Dynamic  in  space  and  time 

•  More  flexible  design 

•  Don't  need  as  many  isopycnals  as 
Z- levels  in  deep  water 

•  Developed  by  NOPP  HYCOM 
consortium 


Can  HYCOM  do  both  global  and  coastal  domains? 


We  will  nest  coastal  HYCOM  and  NCOM  within  regional  and/or  global 
HYCOM  with  progressively  finer  horizontal  and  vertical  resolution 


•  Sub-kilometer  horizontal  resolution  is  the  goal 

•  Optimal  vertical  coordinate  over  shelf  and  shelf  slope 
(isopycnals  over  the  shelf  when  water  is  stratified?) 

•  Sensitivity  to  nesting  ratios,  number  of  nests,  type  of  boundary  condition 

•  Dynamical  impact  of  increased  horizontal  and  vertical  resolution 

•  On-line  nesting  capability  for  multiple  coastal  nests 


Run  on  massively  parallel  supercomputers 


Mississippi 
Bight  domain 


North  Atlantic 
(NATL)  Domain 


Intra-Americas 
Sea  (IAS) 
Domain 


Mississippi  Bight 


Test  cases  are  for  the  Mississippi  Bight 

•  Primary  focus  area  for  two  NRL  projects  using  NCOM  and 
HYCOM  (CO-NESTS)  and  HYCOM  and  ADCIRC  (SEED) 

Large  Mississippi  Bight  domain  using  HYCOM 

•  Nested  inside  an  Atlantic  domain  using  HYCOM,  both  at 
0.08  degree  resolution 

•  Using  off-line,  file-based,  nesting 

•  Makes  the  HYCOM  domain  practical 
(3  processors  instead  of  200+) 

Smaller  Mississippi  Bight  domain  using  NCOM 

•  Always  has  40  fixed  levels  in  the  vertical 

•  Nested  inside  HYCOM  using  archive  files  for  boundary 
exchange 

Two  test  cases  already  performed  as  part  of  CHSSI  project  (Wallcraft) 
•0.08  degree  40-level  HYCOM,  0.08  degree  40-level  NCOM 

•  0.08  degree  26-layer  HYCOM,  0.08  degree  40-level  NCOM 
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Initially  a  HYCOM 
subregion  will  be  set 
up 


hycom  batnymetry 
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26-layer  HYCOM 


Remapped  to  40  levels 


temperature  merid.sec.  87.68w  date:  nov  30,  2002  [03. 2H] 


temperature  merid.sec.  87.68w  date:  nov  30, 2002  [03.2 H] 


salinity  merid.sec.  87. 68w  date:  nov  30,  2002  [03. 2H] 


PLM  remapping  from  26  layer  HYCOM  to  40  level  sigma-z 


Nest  of  40  level  HYCOM  to  40  level  NCOM 


40-level  HYCOM  40-level  NCOM 


temperature  merid.sec.  87.68w  date:  nov  30, 2002  [04. 3H] 


temperature  merid.sec.  87.68w  date:  nov  30,  2002  [94. 3H] 


salinity  merid.sec.  87.68w  date:  nov  30,2002  [94. 3H] 


6  days  after  restart 

Generalized  vertical  remapping  is  the  goal 


THE  END 
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